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We calculate the full I—V characteristics at vanishing temperature in the self-dual interacting 
resonant level model in two ways. The first uses careful time dependent DMRG with large number 
of states per block and a representation of the reservoirs as leads subjected to a chemical potential. 
The other is based on integrability in the continuum limit, and generalizes early work of Fendley 
Ludwig Saleur on the boundary sine-Gordon model. The two approaches are in excellent agreement, 
and uncover among other things a power law decay of the current at large voltages when U > 0. 
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Experimental investigation of transport phenomena in 
quantum impurities is a rapidly expanding field. Typi- 
cally, tiny structures that behave quantum mechanically 
- i.e. , whose level spacing is much larger than all other 
relevant energy scales - are connected to metallic leads 
and a voltage Van across it forces a current to flow. The 
non equilibrium regime that ensues could be achieved 
and measured for structures realized as quantum dots in 
a 2D electron gas [l| or single molecules Q. 

Unfortunately, the theoretical description of these sys- 
tems remains, in spite of numerous efforts, somewhat less 
advanced. One major obstacle is the lack of efficient the- 
oretical approaches to treat non-equilibrium situations 
in the presence of strong interactions. On the analyti- 
cal side, perturbative (Keldysh) techniques are extremely 
difficult to carry out to high orders, while self-consistent 
approximations are difficult to control in the strongly non 
linear regimes. Since impurity problems in the scaling 
limit can be reformulated as ID boundary field theories, 
it is tantalizing to try to use the power of integrability. 
This was done first in Q where full I—V characteristics 
were calculated for the problem of edge-state tunneling in 
the fractional quantum Hall effect. More recently, Mehta 
and Andrei [H have proposed a seemingly different ap- 
proach dubbed the open Bethe ansatz. They also cast 
doubt about the results in Q and their way of coupling 
the model to the reservoirs. While the approach in [J] is 
a priori quite general, it has not, so far, led to quantita- 
tive predictions in the scaling limit because of technical 
difficulties. Besides, since the cut-off in Q is imposed via 
the Bethe ansatz solution, the results cannot be directly 
compared with those of lattice simulations. 

In view of this confusing situation, it would be nat- 
ural to turn to numerical approaches. This is however 
just as challenging. In simulating quantum transport 
one faces the problem that the stationary Schrodinger 
equation is replaced by the time dependent Schrodinger 
equation and therefore an eigenvalue problem is replaced 
by a boundary problem in time. An important conse- 



quence is that in simulations in time (frequency) domain 
one has to send first system size to infinity and then time 
to infinity (level broadening/frequency to zero) . Current 
numerics typically falls into two classes. Either one takes 
the limit of switching on interaction last and starts from a 
noninteracting description, where one can send the sys- 
tem size analytically to infinity. Or one is using non- 
pertubative methods on a finite lattice where the limit 
of time/frequency is interchanged with the limit of sys- 
tem size. Therefore numerical simulations for transport 
properties of strongly interacting quantum system are in 
general prone to conceptual considerations. 

We report in this letter a double step forward. We 
on the one hand extend the approach of j^j to a special 
point in the interacting resonant level model (IRLM), 
where we are able to determine the full I—V characteris- 
tics at vanishing temperature T. We on the other hand 
carry out time dependent DMRG (td-DMRG) calcula- 
tions, and fully confirm our Bethe ansatz predictions in 
the scaling limit. The remarkable agreement between 
analytical and numerical approaches strongly validates 
both, hence dispelling doubts about the issues of reser- 
voir coupling in [3| as well as the feasibility of time de- 
pendent DMRG. On top of this, the physics unraveled by 
our calculations is highly non perturbative and counter 
intuitive, and exhibits, among other things, regimes of 
negative differential conductance. We thus expect the 
IRLM model to become a benchmark for other methods 
in the field of interacting out of equilibrium transport. 

We start with the numerical approach. We apply the 
td-DMRG method [H H, 0, S, ^ to integrate an initial 
state of the nanostructure Hi attached to a left and a 
right noninteracting tight-binding lead 7Yl 
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where t = 1 is the hopping amplitude in the leads, M the 
number of lattice sites, U the interaction on the contact 
link and e c i an on-site gate potential which is set to zero 
(i.e. at resonance). 

At time r = we include a voltage drop by applying 
a potential V SD /2 (-V SD /2) on the left (right) lead which 
we smoothly send to zero on a scale of three sites left 
and right the impurity. We then switch off the voltage in 
the Hamiltonian and time evolve using e in ( WL+Wl ' T (dif- 
ferent procedures to reach the stationary state could be 
considered. We will report on this elsewhere - see also 
[lol]). Since the leads act as a bath the simulation has 
to be stopped after the transit time 7t = Aead/^F, where 
Liead is the leads' length and v F = 2t is the Fermi velocity 
of the noninteracting tight binding leads, since for times 
larger than r t one measures the influence of the bound- 
aries and not the steady state. In addition, since we are 
working with finite leads, our system has a finite size gap 
which leads to finite size induced cos(V^ D r) oscillations. 
This is similar to the oscillations in a Josphson junction 
which are induced by the superconducting gap. By care- 
fully checking for finite size effects and ensuring to take 
enough states per block in the DMRG procedure we can 
extract the current corresponding to infinite leads, the 
details are explained in llj, ll2| . 

In Figure [T] we plot the current / vs. applied source 
drain voltage for a system with t' = 0.5 and a = 0, 
i.e. on resonance. Most of the data was calculated using 
a 96 site system and at least 2000 states per block. For 
comparison we include a reference calculation for U = 1.0 
using 120 lattice sites and 3000 states per block. We 
applied typically 20 to 25 full td-DMRG steps with a 
time step of A T = 0.4 as described in 0. We then switch 
to an adaptive time evolution scheme as described in 0, 
[lil ]. We would like to remark that we perform the time 
evolution using a full Arnoldi type matrix exponential 
during the full and adaptive td-DMRG sweeps without 
any Trotter like approximation schemes. 

Since free fermions provide a non-trivial test for (real 
space) DMRG the U — 0.0 results show that the proce- 
dure is well defined and gives accurate results even for 
large voltages. Switching on the interaction one observes 
in the not too large V SI) regime a broadening of the differ- 
ential conductance while for large interaction the broad- 
ening is absent and one obtains a shrinking of the reso- 
nance width. This is similar to the linear conductance 
vs. gate potential as described in [HI]. For large volt- 
age a negative differential conductance regime appears 
which is maximal at U ~ 2.0 and disappears again for 
U — > 00. We note that DMRG has a tendency to un- 
derestimate the current in the large voltage regime if not 
enough states are kept, especially in the adaptive scheme. 

In Figure [2] we plot the IV curves for t' = 0.2,0.3,0.4 
where we rescaled the / and V^ SD axis in such a way that 
we fit the data to our analytical result using a scale Tb- 

The points nicely sit on a single curve, which shows 
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FIG. 1: (Color online) Current vs. source drain voltage V SD 
for a hybridization of t' — 0.5, and interaction values of U = 
-1.0,0.0,0.3,1.0,5.0, and 10.0. The subscript 'V ('2') of I 
refers to data extrapolated from link to the left (right) lead. 
The lines are guides to the eye, except for U = 0.0 where we 
plotted the exact result for infinite leads. Calculation have 
been performed with 96 sites, 48 fermions keeping at least 
2000 states per DMRG block. In addition data for U = 1.0, 
M — 120 sites and Ncut = 3000 states per block is shown. 



that we are reasonably within the scaling limit. A numer- 
ical fit of the power law decay for t' — 0.2 and U = 2.0 
gives an exponent of 0.47. In the inset of Figure [2] we 
show that the numerically obtained Tb scales are given 
by a Tb ~ {t') i/3 power law. 
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FIG. 2: (Color online) Comparison between analytical and 
DMRG results at the self dual point. The numerical data have 



been fitted using a single parameter Tb = c(t 
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Exponents for different values of U are represented on 
figure [3l For all positive values of U we find a power law 
decay of the current at large voltage, except for t' = 0.5 
and U = 0.3 where the crossover scale is beyond our volt- 
age regime. We also find that the exponent does not vary 
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monotonously with U, and reaches its maximum around 
U = 2 where it is very close to 1/2, and approaches 1/2 
when lowering t' , i.e. in the scaling limit. 




FIG. 3: (Color online) Exponents in the negative differential 
conductance regime for t' = 0.2, 0.3, 0.4, and 0.5. 



We now turn to a study of the model in the con- 
tinuum limit. It is convenient to first unfold the 
left and right leads, obtaining in this way Hl = 
— i J2 a =i 2 I dx ipl i d x ip a a free Hamiltonian describing 
two infinite right moving Fermi wires, and 

Hi =^Oi(0) + V 2 (0)) rf+h.c. (3) 
+ U C (: ^|^(0) : + : ^ 2 (0) : ) (cftd - |) + e d Sd 

In the following, we focus again on the resonant case, 
i.e. ed — 0. Note that in ([3]) the Coulomb interaction 
involves normal ordered charge density on the dot, as in 
([2|) . The subscript " c" indicates that the precise relation- 
ship between the lattice and continuum limit parameters 
depends on the regularization scheme. 

We switch to the language of the anisotropic Kondo 
model by representing the impurity degree of freedom by 
a spin one half, <$ = rjS + , S z — Sd — |, with rj a Majo- 
rana fcrmion. The fcrmions in the wires are bosonized as 
ib a = -4= e Jv/ ^ ¥ ' a . We then introduce new bosonic fields 

ra V2tt _ _ 

(fc = (if i +(^2)/v / 2 and cj) — (fi — f2)/V2, that represent 
the total and relative charge in the wires respectively. 
The unitary transform IA — exp (ia^/2S z if c (0)) is then 
applied. Choosing a = U c /^/tt to cancel the boundary 
interaction along S z , we arrive at: 



H = H (ip c ) + H ((f>) + Hi 
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where Tto(X) = dx (d x X) 2 is the Hamiltonian of a 
free (chiral) boson, and the piece involving the impu- 
rity degrees of freedom is Hi = ^ [V(0)O(0)S+ + h.c], 

with y = e ^<=, O^K 1 e tV ^^ + K 2 e- lV ^ ,t ', and K a = VVa- 
The exponent of the vertex operator is (3 C = \/2~%(l — ^ SL ). 
We note that there is a special value of the coulom- 
bic repulsion (U c — n in our scheme) where a remark- 
able simplification occurs: this exponent vanishes, and 
the field tp c decouples from the impurity. At this point, 
and this provides a universal characterization thereof, 
the anomalous dimension of the boundary perturbation 

D — t + (| — j 51 ) reaches its minimum value D = \. At 



Although we do not have full analytical solutions for 
general values of U c , some qualitative results can easily 
be obtained. First, we expect that, provided U c > 0, 
the current vanishes at large voltages like a power law 
Iq(V) ~ 5=1 — 2D, D the scaling dimension of the 
perturbation (this conclusion was reached in discussions 
with B. Doyon [HI). From the foregoing discussion this 
leads to an exponent b = (2 — ^) We emphasize 
that this only holds in our regularization scheme. The 
numerical values of U c where b reaches its maximum and 
where b vanishes do not have to be the ones observed 
in the lattice model; in particular U c = 2ir for us cor- 
responds presumably to U = oo where it can be argued 
that the model is equivalent to U — 0. Nevertheless, it 
is usually expected that different regularizations do not 
change the qualitative nature of the results. The pre- 
diction for the exponent is indeed in full agreement with 
numerics and the curve on figure [3l We see in particu- 
lar that the exponent b does exhibit maximum value at 
b = h where D = i by scaling. This allows us to identify 
the self dual point in the continuum limit with the value 
U ~ 2 in the lattice model. We shall now see that this 
self dual point is also amenable by the Bethe ansatz. 

The IRLM in equilibrium is solvable: the correspond- 
ing basis of the Hilbert space provides quasiparticle ex- 
citations which scatter diagonally (ie, without particle 
production) across the impurity. In general however, the 
operator enforcing the voltage drop across the impurity 
^ J dx(ip\ipi — i\)\i\)2) = ^f-Q is not diagonal in the 
quasiparticle basis. This makes the construction of exact 
scattering states along the lines of 0] a seemingly impos- 
sible task. At the special value of the coupling U c = tt 
however, things simplify drastically. The hamiltonian 
can be mapped onto the boundary sine-Gordon theory 
for the field <fi after folding across the impurity (it can be 
shown that the cocycles can be discarded at this stage). 
This theory is known to be integrable, and leads to a de- 
scription in terms of solitions, antisolitons, and two kinds 

of breathers. The charge Q= \ j d x (j) is not conserved 
by the interaction, but acts diagonally on the quasipar- 
ticle basis. Following the arguments in 0] this allows for 
an exact calculation of the non equilibrium current at all 
values of T and voltage for the IRLM at XJ c =ix. We will 
simply give here the result at T = 0. 

In this case, the in-state involves one type of quasi- 
particle - for Vsd > say it is the antisolitons. We 
parametrize their energy (equal to the momentum) as 
p = wnere m is an arbitrary mass scale, and A 

the rapidity. The antisoliton filling fraction thus reads 
f-(\) = Q(A - X) (9 is the Heaviside function), A be- 
ing a rapidity cut-off. Defining p = nf- , where n is the 
density of allowed states per unit of length and rapidity, 
the current reads (we use units in which e = h = 1) 



this point, the model exhibits a certain self-duality [14 1 
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Here, Xb is a rapidity encoding the crossover energy scale 
in the problem T B = r fe Xs . From scaling, T B oc {t' c ) 4/3 . 
Like for the Kondo temperature, different ways of defin- 
ing Tb are possible. We follow here the definition given in 
[10 ] which is related to universal coefficients in the low-T 
expansion of the linear conductance; we will not need it 
in the following anyway. The density n - or equivalently 
p - follows from non trivial Bethe ansatz quantization 
rules which involve the antisoliton-antisoliton scattering 
matrix. The Wiener Hopf technique allows to obtain 
the Fourier transform p(w) = J_ 
form. One finds (see Ref. [171]): 



dX e lulX p{\) in closed 
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e iujA and e z 
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The 



With G(uj) = V°7< r ^to\ r /l+j M T 

relation between the cutoff A and the applied voltage can 
be expressed through the Fermi momentum p F — ^e A — 

r(2/3) ' P° wer expanding the denominator of the 
integrand in ([5]) yields an explicit series representation 



with natural expansion variable V SD — 4^/^(2/3) tS" 
Depending on whether V SD is small or large one has: 
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For comparison, recall that in the free fermion case 
(U=U c =0) one has I (V SD ) = ?f Atan(^). We see in 



(|7l8p that Iq/Tb is a function of V SD only, so that match- 
ing the numerical results in the scaling limit requires a 
single common rescaling of I and Vsd- The T = current 
is depicted in figure [2] (solid line) . It decays as a power- 

1 . 1 u T 3 3/4 r(2/3) 9/2 T1 3/2 T ,-l/2 

law at large voltage, J ~ %r TwV" r s ^ SD ■ 

We can understand more precisely the origin of the 
large V SD decrease of the current by looking at the non- 
equilibrium density of states for anti-solitons: p(p) = 

L'lfi^ wHln T? (n-\ - 3 £(l) (~T 



F p (^)6(p F -p) with F p (a: 



87? 1 Sn>0 (2n+l)! 

At small momentum p <C 



r(^)r(M)(^)^ 

VIjd - where the resonance (of width Tb) forms - the den- 
sity of states vanishes as a power law, p(p) ~ y/p/V SD : 
this depletion of the sea close to zero energy results in an 
extinction of the current. The situation is clearly con- 
trasted to what happens in a free theory (U — 0), where 
the density of states is constant provided p < p F — V SD . 

We note that the formula b(U c ) suggests a power law 
divergence of the current at large voltage in the scaling 
regime for U c < 0. This seems confirmed by our numerics. 

In conclusion our work provides what may be the first 
example of transport properties in an interacting one di- 



mensional system out of equilibrium that can be calcu- 
lated both analytically and numerically, with excellent 
agreement between the two approaches. This should pro- 
vide a most useful benchmark for the variety of other 
approaches being currently proposed. Our results ex- 
hibit remarkable physics: the negative differential con- 
ductance at large voltage seems a truly non perturbative 
behaviour, with unclear physical origin: a possible expla- 
nation could be that once we are in the tail of the conduc- 
tance curve a voltage drop at the impurity builds up since 
differential conductance is now smaller than unity. This 
voltage drop may destroy the interaction based renormal- 
ization of the conductance enhancement. 

The IRLM provides a perfect laboratory to explore 
other challenging questions, such as the DC shot noise, 
or the effect of coupling of the baths on the stationary 
properties. We hope to get back to these soon. 
Acknowledgments: we are grateful to N. Andrei, B. 
Doyon, F. Essler and P. Mehta for numerous discussions 
and encouragments. HS was supported by the ESF net- 
work INSTANS. The DMRG calculations have been per- 
formed on HP XC4000 at Steinbuch Center for Comput- 
ing (SCC) Karlsruhe under project RT-DMRG. 



[9] 
[10] 

[11] 
[12] 



[13] 

[14] 
[15] 
[16] 
[17] 



See for example D. Goldhaber-Gordon et at, Nature 391, 
156 (1998); D. Goldhaber-Gordon et al, Phys. Rev. Lett. 
81, 5225 (1998). 

W. Liang et al, Nature 417, 725 (2002). 

P. Fendley, A.W.W. Ludwig and H. Saleur, Phys. Rev. 

Lett 74, 3005 (1995). 

P. Mehta and N Andrei, Phys. R ev. Lett. 96, 216802 
(2006). Erratum [cond-mat /0703426| 
S.R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. 
B 48, 10345 (1993). 

S.R. White and A.E. Feiguin, Phys. Rev. Lett. 93, 076401 
(2004) . 

A. J. Daley, C. Kollath, U. Schollwock, and G. Vidal, J. 
Stat. Mech.: Theor. Exp. (2004) P04005. 

Density Matrix Renormalization - A New Numerical 
Method in Physics, ed. by I. Peschel, X. Wang, M.Kaulke, 
and K. Hallberg, Springer (1999); R.M. Noack and 
S.R. Manmana, AIP Conf. Proc. 789 (2005) 93-163. 
P. Schmitteckert, Phys. Rev. B 70, 121302(R) (2004). 

B. Doyon and N. Andrei, Phys. Rev. B 73, 245326 (2006). 
G. Schneider and P. Schmitteckert, condmat/0601389. 
G. Schneider and P. Schmitteckert, in High Performance 
computing in Science and Engineering '06, edited by 
W.E. Nagel, W. Jager, and M. Resch (Springer Verlag, 
Berlin, 2007) pp. 113-126. 

D. Bohr and P. Schmitteckert, Phys. Rev. B 75, 
241103(R) (2007). 

A. Schiller and N. Andrei, cond-mat/0710.0249 

B. Doyon, Phys. Rev. Lett. 99, 076806 (2007). 

E. Boulat and H. Saleur, Phys. Rev. B 77, 033409 (2008). 
P. Fendley, A.W.W. Ludwig and H. Saleur, Phys. Rev. 
B 52, 8934 (1995). 



